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ON THE DEVELOPMENT OF THE PEETUEBATIVE FUNCTION 
IN THE THEOEY OE PLANETAEY MOTION. 

By E. T. A. Innes. 


(Eead April 19, 1916.) 


In the Transactions , yoI. ii, part 3, pp. 301-317, 1911, I have a paper 
upon the Newcomb operators used in the algebraical development of the 
elliptic perturbative function. The present paper deals with a further 
extension of the uses of these Newcomb operators. At best the develop¬ 
ment of the perturbative function is a very cumbersome business, and 
reference must be made to Newcomb's original paper in vol. v. of the 
Astronomical Paj^ers of the American Eyhemeris for a full statement. On 
p. 15 Newcomb quotes his remarkable theorem by which the operators, 
acting on the development in powers of e, combine with the operators 
acting on the development in powers of e lf e and e 1 being the eccentricities 
of the inner and outer planets respectively. The operators connected with 

e he writes as n”, those with e l as nj"/, and he proves that 

u 7j' — u j x n 0 , r 

These operators consist of a string of differential coefficients of the Lap- 
lacian functions hi or a'A^ 

Thus we have— 


8Ili’o = [8+ (4i* - 2i)D +(-2i- 1)D 2 - D"]6 ; , 

in which— 

r> d 

da 

Operators which are factored by powers of a 2 = sine squared half mutual 

k 

inclination of orbits are indicated by the affix Jc thus : {\ n > n> . 

Unfortunately Newcomb only shows how the operators can be combined 
symbolically, and therefore in numerical work it is necessary to consider 
the eccentricities of both orbits simultaneously. 
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It appears, however, that an extension can be made. It will be proved 
that if we possess a development of the disturbing function in which one 
of the eccentricities is supposed to be vanishingly small (in fact zero, but 
with a place of perihelion), we can include the effect of the second eccen¬ 
tricity by a series of simple operations. Thus, starting with a development 
with e = 0 and e x = e l we can arrive at the complete development for 
e and e v 

Newcomb writes any term of the development as follows : 


eie/a- k Vj, f (i)C(i'(j' - ig + &),* 

in which— 


k k k 

P/,/ = Ilf, + + e\Uj’J, + terms factored by e 4 , e l e\, e 4 , etc. 

In applying the new method it is a matter of choice which eccentricity is 
assumed to be zero. If each in turn is taken as zero and the two partial 
developments completed, then an independent check is furnished. Generally 
as flo,/ is larger than Uj it will be better to assume that e, the eccentricity 
of the inner planet, is zero. 

Implicitly the method not only computes the development, but also its 
derivatives; but as these are actually wanted, it is so much to the good. 

If we suppose that e 1 is zero, we will indicate the fact by replacing the 
upright P used by Newcomb by a slanting P; if c is zero, we must add a 
dash so as to distinguish between— 

Po,o = no,o + + etc., 

and— 


Po,o = bIo,o + d- e fblo’o d" e t c * 

Assuming now that is zero, any term in the development takes 
this form— 
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k 



We cannot write down the general expression for Pj, but we know it 
contains a series of terms— 

k k k 


n 7 . 4- e*n 


2,0 

j, 0 


+ e^r/o + etc. 


* Leaded C indicates “ cosine/' 
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The question is if we operate on Pj with say the operator no.i = 
1 .o 

o(2 i + 1 + D), which reduces ]J.f,o to etc., so as to obtain P/,i, can we 

z . k 

we include P±j_ 2 fc in this operation? We can see that this is so for 

k o 

j positive because the operator 17,(1) is equal to the operator n,(i) if i is 
replaced by i—lc. This will also include the negative suffices, because we 

k k 

pass from fl^ to 17 by replacing i by — (i — 21c). 

Hence— * * 

P - 7 — 2 »( — i — &) — P; + 2 n (f — h). 

Similarly for the external planet the operator n 0 /W i s equal to 

17 (i + 2&), 

and— k k 

Po,-y- 2 »( — i — lc) — Po,y+ 2 »(f — &). 

Hence when we have reduced each term to the form— 

pc& coy + ig) — ps& soy + ig)* 

PCQ' C(i'g + ig) — PSfl'SOY + 0) 

2 2 
PCd"C0y+ ig) — PS6"S(iy + i</) 

we can operate on the sum so as to include the effect of the eccentricity 
of the outer planet—and vice versa. 

As already mentioned, an absolute check is obtained by considering 
firstly one eccentricity nil and then operating, and then considering the 
other nil in its turn and operating. The resulting sums for each term must 
be identical, and they will be so if all sensible terms are included and the 
calculation is correct. 

Each operation introducing e removes the term from the argument 
ig -I- ig to the next argument i r g' + (i -f \)g, and the introduction of 
e n removes it n times. 

As an example we give a portion of the process as applied to some of 
the terms which go to make up the great inequality of Jupiter and Saturn. 

Suppose, firstly, that Jupiter’s orbit is treated as being without 
eccentricity, 

We have then — 

e i-Po ,3 (2)C[5 g’ — 2 g + 2 (w’ — w ) +0(w' + «>)] 

4- ey-Po.i (3)C[5flr' — 2g 3 O' — w) + 1 O' + w)~\ 

4- c^HVi(4)C[%' - 2 g + 4<>' -w) + 20' + w)] 

-f terms of the 9th order. 


* Leaded S indicates " sine.” 
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Into these coefficients, and until they are insensible, we introduce the 
powers of e by operating firstly with H 2 = ±( — 16 4 - D + D 3 ), then with 
n<!> and so on. ‘ 

If these results are added together, the effect is to change P into P. 

The next coefficient will arise from the term— 

e?Po, 2 (B)C[ 5 ^f / — 3g + 3 Qw' — w) + 0(w' + ui)] 
i 

+ o- 2 Po.o(4)C[5^ / — 3 g + 4 (w' — w) 4 - 1 (w' -f w)] 

4 - etc. 


In this case the succession of operators is-ellJ + e*Y\\ -r e 5 Ili+ etc., and 
so on. 

Similarly the term starting with— 

Po, 4 ( 1 )C[%' — g + O'— w)], 

when operated on with ell 1 -! + e 3 n^i + etc., 
will give the coefficient P_i. 4 (l) etc., etc. 

In practice each term PC [5/ — ig + 0 ] will be brought to the form— 

Po,nOC[5/- ig ] - P 0jn S6S[5g' - ig ] 
i 

and all the coefficients PC0, PC0', etc., summed so that we have finally— 
2PoyC0C[5/ - ig -]SP Of *'S0S[5/ - ig-]. 

Then the operator which starts with e n U^ + e w + 2 n*J 2,0 + etc., acting 
on this sum will yield us— 


SP, f X0C[5 g' - (t - n)g] - SP W# „/S0S[5/ - (i - n)g ], 
which is complete in P M w >. 

In thus making a complete expansion of the perturbing function and its 

derivatives, the saving of labour will be great. In Newcomb’s and Le 

Terrier’s expansions each coefficient requires a separate and sometimes very 

lengthy computation, whereas, as just demonstrated, blocks of coefficients 

are computed in one process. The machinery by which the other eccentricity 

is introduced is uniform throughout and simple. Moreover, it is more 

inclusive. In algebraical expansions of the perturbing function, the com- 

putor must decide at what limit he will stop. In Le Verrier’s theory of 

Jupiter and Saturn, he adopted roughly the 4th and 5th powers of the 

eccentricities. In other words, he was satisfied with such an expansion as_ 

^n 0 0 ;l + e^Uoii + e ~ e ifi 2 ;J + 

The new method starting with— 

P n 0 ’ 1 4- P 3 n 0 ’ 3 4- ^n°> 5 
e i u o,i i e i 1 J o,i t- V J o,i> 

which is the same order of approximation, will at the first operation 
introduce— 
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The next operation indicated by flo o introduces— 
e i e 1 Uo’} + 


4,5 
‘ V 


and so on, and of these Le Yerrier only included the first of each. The 
point is, that with exactly the same labour as Le Yerrier used for the lower 
terms, we introduce in addition a crowd of higher terms, some of which may 
be sensible. 


Newcomb, on pp. 27-32, gives the expressions for the operators n* an d 
n q y so far as regards development to the 8th powers of the eccentricities— 
these or the modification of them already given in Transactions , vols. ii, 
part 3, and iii, part 2,1911 and 1913, must still be used; but the tables for 
II};}/, which Newcomb gives on pp. 32 to 42 are not required. Similarly for 


the II. The number of pages involved shows pretty closely the ratio of the 
saving—viz., as 4J pp. are to 10 pp. 

Johannesburg, March 27, 1916. 


P.S. —The process outlined above can be generalised. If we have any 
development involving the distances of two planets from the sun, and 
suppose both eccentricities to be zero, then the effects of the eccentricities 
can be introduced bit by bit so far as they are sensible. Thus, if with 
eccentricities zero, we have a term 

/(a)COY - ig + 0) 

then the Newcomb operator Pjj< will introduce the effect of ej e x j' and 
change the angle from— 

i'g' - ig 4- 0 

to— 

O' + /) g' - 0 - i) g + 9- 

A numerical example of this new view of the use of the Newcomb operators 
will be found in the Astronomical Journal for June 30, 1916. 

September 9, 1916. 



